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Abstract 

We consider quantum field theory on a spacetime representing the Big 
Crunch/Big Bang transition postulated in ekpyrotic or cyclic cosmologies. 
We show via several independent methods that an essentially unique match- 
ing rule holds connecting the incoming state, in which a single extra dimension 
shrinks to zero, to the outgoing state in which it re-expands at the same rate. 
For free fields in our construction there is no particle production from the 
incoming adiabatic vacuum. When interactions are included the particle pro- 
duction for fixed external momentum is finite at tree level. We discuss a 
formal correspondence between our construction and quantum field theory on 
de Sitter spacetime. 
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I. INTRODUCTION 



Despite its overwhelming phenomenological success, the standard big bang cosmology 
is clearly incomplete. Its gaps and paradoxes provide some of the most powerful clues to 
fundamental theory that we possess. Indeed, it is increasingly evident that the real measure 
of success for string theory and M theory will be how well they face up to the challenges posed 
by cosmology. Perhaps the greatest challenge is that of describing the initial singularity, a 
moment of infinite density and curvature occurring some fifteen billion years ago in our past, 
a basic puzzle not resolved by cosmic inflation. 

The initial singularity is often associated with the problem of the 'beginning of time'. 
But the only thing one can legitimately infer from the existence of the singularity is that 
general relativity is incomplete. Rather than have time 'begin', which is a truly paradoxical 
notion, or to work with imaginary time formulations, it seems reasonable to explore the 
alternative possibility that time may be continued back through the singularity, and even 
arbitrarily far into the past. Such a view is consistent with what is known so far in string 
and M theory. Spatial geometry and topology are only approximate concepts, as evidenced 
by orbifold backgrounds,! and allowed topology changing processes.! However, time is built 
in, in a fundamental role, and there is no evidence so far that it is allowed to 'begin' or 'end'. 

Recent attempts to construct cosmological scenarios employing 'brane world' construc- 
tions from M theory, and string theory have led to a re-examination of these issues. The 
'ekpyrotic' scenario^ in which a brane collision is supposed to be the origin of the hot big 
bang, and its 'cyclic' version^ in which such collisions occur periodically into the infinite past 
and future, provide alternate approaches to the classic cosmological puzzles conventionally 
addressed by inflation. In the cyclic model, the flatness, homogeneity and isotropy of to- 
day's Universe is explained as a consequence of an epoch of vacuum energy domination in 
the previous cycle. And the density perturbations needed to seed structure formation were 
generated by an inter-brane attractive force near the end of the last cycle. An important 
precursor of these ideas was the 'pre-big bang' model of Veneziano et alB 

The ekpyrotic and cyclic models rest for the most part on conventional low energy ef- 
fective field theory and gravity. One key event cannot be described within that approach, 
namely a collision between the two end-of-the- world boundary branes (or 'orbifold planes'). 
In the four dimensional effective description this event appears to be unavoidable, since the 
four dimensional effective scale factor is initially contracting. The four dimensional fields 
appearing in the theory have positive (and growing) kinetic energy and this means, through 
the Friedman equation that the contraction cannot be reversed. Within a finite time one 
reaches a 'big crunch' singularity dominated by scalar kinetic energy, an event which appears 
at first sight to be irredeemably singular. From the the higher dimensional viewpoint the 
situation is more optimistic. The geometries of the branes are regular at the collision and 
the density of matter on the branes is finite. The five dimensional Riemannian curvature is 
finite everywhere away from the singular point. In fact, the only sense in which the higher 
dimensional geometry is singular is that the fifth dimension shrinks away to zero size.! 
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FIGURES 




FIG. 1. The compactified Milne universe. On the left is two dimensional Minkowski space. 
The Lorentz invariant coordinate t satisfying t 2 = T 2 — Y 2 is constant on the dashed surfaces, 
which provide a spacelike foliation of the causal future and past of the origin. These surfaces 
are parameterized by a coordinate y. Identifying y with y + L compactifies space to produce the 
spacetime on the right, consisting of two Lorentzian cones joined tip-to-tip at t = 0. If the circular 
sections of these cones are orbifolded by a Z2, then the two fixed points of the Z2 are two tensionless 
branes which collide and pass through one another at t = 0. 

It is crucial for the cyclic scenario, as currently formulated, that a satisfactory method 
be found for passing through the singularity corresponding to the collapse of the extra 
dimension. In particular, the issue of matching the density perturbations across this sin- 
gularity has been a matter of fierce debate.0 A matching rule was proposed in Ref. || 
according to which the growing mode scale invariant density perturbations developed in the 
pre-collapse phase are transmitted across the singularity. But it is also possibleQ to match 
in such a way that only the decaying mode is present in the final state. Interesting papers 
have subsequently appeared suggesting geometrical methods of regularizing the singularity,! 
or employing scalar fields with a negative kinetic term to do so.M However, none of these 
methods yet yields a completely unambiguous result for the case of interest in the ekpyrotic 
or cyclic scenarios. We hope that the method developed here on more fundamental grounds, 
when extended to include gravitational backreaction, will be applicable to the cosmological 
case. 

Ultimately this issue must be dealt with by string or M theory. Indeed, regardless of the 
ekpyrotic or cyclic scenarios, there are good reasons for believing that this type of singularity 
must be resolved if string theory is to make sense. The shrinking of the extra dimension can 
be accurately described using a slow motion (moduli space) approximation, which remains 
valid all the way to zero size. The low energy moduli of string theory and M theory are 
believed to be fundamental, more so even than the actions and Lagrangians they are derived 
from.0 For example, these moduli are the parameters which interpolate from one corner of 
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M theory to another. The shrinking away of one extra dimension, in finite time, seems to be 
perfectly allowed in string theory, either if it is one of the nine string theory dimensions, 
or if it is the tenth spatial dimension associated with M theory.! In the former case, the 
string coupling is constant and may be taken to be arbitrarily small, so stringy interactions 
should be negligible. In the latter case, the string coupling vanishes as the extra dimension 
shrinks away. Non-perturbative effects should, in this case, vanish even more rapidly than 
perturbative effects. Thus it is hard to see what would prevent this process. The question 
which must then be faced is: What happens next? 

The moduli space description, and the higher dimensional picture, both lead to a natural 
continuation,! illustrated in Figure 1. The extra dimension contracts to zero at a certain rate 
but immediately reappears at the same rate. In the brane picture, the two branes collide 
and pass through one another, a behaviour familiar from BPS solitons in other contexts. 
If the collision occurs at finite speed, one expects some associated particle production and 
consequent back-reaction. 

In this paper, we take modest steps towards our eventual goal of a calculation of the 
consequences of a collision between boundary branes in M theory. There are significant 
technicalities to be faced even at the level of quantum fields, which is all that we shall 
discuss here. We shall propose a method of obtaining a unitary quantum field theory on 
the spacetime illustrated in Figure 1. Within free field theory, in our construction there 
is no particle production in passing from the big crunch to the big bang phase. However, 
once interactions are included, particle production occurs. For fixed external momenta, the 
particle production at the big crunch/big bang transition which is well defined and finite at 
tree level. It exhibits a power-law fall-off at high momenta which we argue would likely be 
replaced by exponential fall-off in string theory. 



II. MILNE AND COMPACTIFIED MILNE 

The spacetime we are interested in is a subspace of d + 1 dimensional Minkowski space, 
a trivial solution of d + 1 dimensional general relativity or supergravity. The line element is 

ds 2 = -dT 2 + dY 2 + dx 2 , (1) 

where we adopt units in which the speed of light is unity. We shall refer to Y as the fifth 
coordinate, having in mind the picture that three of the d—1 coordinates x should provide 
the spatial dimensions of everyday existence, with the remainder compactified for example 
on a torus or orbifold. For example in eleven dimensional M theory, d = 10 and six of the 
x dimensions would be taken to be compact. 

The line element (fj) may be rewritten in terms of new coordinates defined by T = 
tcosh(H 5 y), Y = tsmh(H 5 y), where — oo < t < oo and — oo < y < oo cover the causal 
future (t > 0) and past (t < 0) of the origin Y = T = 0. We have here introduced the 
parameter H 5 , with dimensions of inverse time. In these coordinates, ([I]) becomes 

ds 2 = -dt 2 + H 2 t 2 dy 2 + dx 2 , -oo < t < oo, -oo < y < oo. (2) 

The space comprising the causal future and past of the origin, and its light cone T = ±Y 
is what we shall define to be the Milne universe M. . The complement of M. in Minkowski 



4 



space comprises the two Rindler wedges |T| < \Y\ to the left and right of the origin in Figure 
1. 

The second step in obtaining the compactified Milne universe Aic is to compactify the y 
coordinate into a circle. Because Ai is invariant under translations, in the quantum theory 
there exists a unitary operator A(H^L) implementing y — > y + L, which is just a boost of the 
original T, Y coordinates on Minkowski space, with rapidity H§L. The coordinate t, which 
is the time in the Milne universe, is invariant under this operation. Let Y(H§L) denote 
the discrete group generated by A(H 5 L). Then we define Aic to be Ai/T(H 5 L), i.e. the 
spacetime 

ds 2 = -dt 2 + H 2 t 2 dy 2 + dx 2 , -oo < t < oo, < y < L, (3) 

where y = and y = L are identified. We see that the parameter H 5 is just the rate of 
expansion or contraction (or 'Hubble constant') of the fifth dimension. The space Aic is 
not a manifold, since it is not Hausdorff at t = 0. But of course this is precisely the point 
of interest to us. 

So far branes have not entered. We may however further reduce the circle < y < L, by 
identifying its upper and lower halves under the Z 2 symmetry y — > L — y. Quantum fields 
may be decomposed into components which are even or odd under this operation. The two 
fixed points of the Z 2 symmetry, y = and y — L/2 can then be viewed as (zero tension) 
orbifold planes, which collide and pass through one another at t — (Figure 1). 

We shall also be interested in studying quantum fields in this background from the point 
of view of the dimensionally reduced rf-dimensional theory. Writing the d + 1 dimensional 
line element as 

dg 2 = e 2^(d-2)/(d-l) dy 2 + e -WVfFWVgWdx' i dar, (4) 

the d + 1 dimensional Einstein action reduces to that for d dimensional gravity with a 
massless, minimally coupled scalar field 0. (We adopt units in which the coefficient of the 
Ricci scalar in the d dimensional Einstein action is ~.) The solution Aic x -R d_1 is now 
re-interpreted as a cosmological solution in which the d- dimensional Einstein-frame metric 
g (d) = a 2 r)^ with scale factor a oc \t\ 1 ^ d ~ 2 \ and <p = J(d - l)/(d - 2)ln|# 5 t|. 

It is clear that gravitational waves travelling in the noncompact directions are minimally 
coupled both in the d+ 1 dimensional description, and in the d dimensional description since 
the powers of <fi in (^j were chosen to obtain Einstein-frame gravity in the reduced theory. 
It is straightforward to check that a scalar field which is minimally coupled in the d + 1 
dimensional theory is also minimally coupled in the d dimensional theory. This means that 
for the background Aic x the dimensionally reduced action for a minimally coupled 

scalar <p is — | / ^J—g^g^^d^ipd^ip = — \ J dt^^d^ipdyip, for any d. 

III. FREE FIELD BEHAVIOUR ON M c 

Let us now describe the behaviour of free fields on Aic x R d ~ l ■ Expanding the fields in 
plane waves e lk ' x on R , modes of momentum k aquire a mass squared of k 2 in their two 
dimensional {t,y) action or equations of motion. The two dimensional line element is just 
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—dt 2 + H 2 t 2 dy 2 , which is conformally flat, with a conformal factor which vanishes at t — 0. 
The two regions t < and £ > of M. are each conformal to an infinite cylinder labelled by 
a conformal time r±, defined by ifst = ±e ±H5T± , in the two cases. The line element in these 
coordinates is then 

ds 2 = e ±2H ^{-dr 2 ± + dy 2 ), (5) 

where — oo < t± < oo on each cylinder. The conformal factor vanishes as t tends to zero. In 
two dimensions the kinetic term for a scalar field is conformally invariant, and hence does 
not see the conformal zero. But a two dimensional mass term vanishes like \t\. Therefore, 
in the limit t — > 0, all field modes behave as those of a massless two dimensional field on an 
infinite cylinder. Modes with nonzero y-momentum k y oscillate an infinite number of times, 
as e ±lkyT± or \t\ %ky / H6 , as r± tends =foo. On the other hand, modes with k y = instead 
evolve linearly in r±, which means that they generically diverge as log|t| as t — > 0. 

The problem of defining a quantum field on Ate is that of matching the modes across 
t — 0, from their asymptotic behaviours as t tends to zero from above or below. Since the 
modes either undergo an infinite number of oscillations, or are logarithmically divergent, 
this matching is quite subtle. Let us discuss the k y = modes in more detail. The general 
solution for the k y = modes behaves as (p ~ A+B \n(\H 5 t\) as t approaches zero, with A and 
B two arbitrary constants. As we approach t — the scalar field diverges logarithmically 
but its canonically conjugate momentum ir = \H5t\1fi tends to a finite value H§B. Our 
problem is then to match a general incoming solution tpo(t) = A~ + B~ \n(\H 5 t\), t < 0, to 
the corresponding solution for t > 0, (fo(t) = A + + B + \n(\H 5 t\). A crude approach would 
be to simply cut the spacetime off at t = ±5, identify the field and its conjugate momentum 
on the two surfaces, and take the limit of small 5. Since the momentum is time-independent 
we obtain B + = —B~, independent of 5 as 5 — > 0. But matching the field yields the cutoff- 
dependent result A + = A~ + 2B~ ln(H 5 S), which implies that for any regular in state, the 
amplitude of the mode functions generically diverge logarithmically with the cutoff. If we 
were to accept this result at face value, the number of particles produced would diverge as 
the square of the logarithm of the cutoff. It is tempting to think that this is a consequence 
of the unphysical sharp cutoff and that a smoother regularization prescription might remove 
the divergence. However, a smoother cutoff, such as replacing \t\ in the action by y/t 2 + S 2 
(geometrically, this amounts to replacing the singular spacetime M.q with an hourglass, 
whose waist has circumference 8H 5 L), leads to exactly the same logarithmic divergence as 
5^0. More sophisticated methods must be sought for making quantum field theory on 
M.c wen defined, as we now explain. 
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FIG. 2. Our first method for constructing quantum fields on Ai, illustrated in a conformal 
diagram of Minkowski space. The unitary map from the past light cone of the origin, t = 0~, to 
the future light cone t = + is defined by free field evolution across the Rindler wedges to the left 
and right of the origin. Using this rule we obtain a unitary theory on the Milne universe M, which 
may then be compactified into the space A4c shown in Figure 1. 

IV. QUANTUM FIELD THEORY ON M c 

We shall describe several different constructions for quantum fields on A^c, which all 
yield an essentially unique result. 

The first method is based on Figure 2. We use the embedding of M.q in Minkowski 
spacetime to define the map from t = 0~ to t = + . This is possible as long as one or more 
of the x directions are noncompact, because in this case, the corresponding momenta k are 
continuous and k = is a set of measure zero. From the two dimensional standpoint, this 
means that all modes are effectively massive. And for massive fields, free field evolution 
provides a unitary map between the past light cone of the origin (t = 0~) and the future 
light cone (t = + ), because no information can be carried off to null infinity J + . The 
coordinate t analytically continues to a spacelike coordinate in the Rindler wedges, and as 
one follows the trajectory plotted in Figure 2, this coordinate runs from zero to a finite value 
then back to zero. So in effect a 'clock' measuring t makes no progress whilst the trajectory 
is outside the region M. of interest. 
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FIG. 3. Integration contours used to define the positive and negative frequency modes on the 
entire Milne universe M. 

This first method may also be viewed as a certain analytic continuation in the complex 
t-plane (Figure 3). The field equation is analytic in the original Minkowski coordinates T 
and Y so the global solution may be obtained unambigously by analytic continuation in 
those coordinates. We shall show that this corresponds to the continuation of the positive 
and negative frequency mode functions and from negative to positive t, illustrated 
in the diagram above. The positive frequency mode functions so defined are analytic in the 
lower half t-plane, and the negative frequency mode functions are analytic in the upper half 
t-plane. The quantum field, being a sum of the positive and negative frequency modes, is 
continued in this mixed fashion across t — 0. This analytic continuation method is actually 
more fundamental than continuation across the Rindler wedges, because it does not involve 
those unphysical regions. This is an important distinction when we introduce interactions. 
There is an ambiguity (for example about what the mass used in the free field propagation 
should be) in the Minkowski space continuation, but no corresponding ambiguity in the 
method illustrated in Figure 3. 

Nevertheless it is interesting to discuss in more detail how the two methods correspond, 
for free fields. The coordinate t continues to a spacelike variable t = ±is in the Rindler 
regions, where the line element is ds 2 — s 2 dy 2 , and y is now timelike. So in the Rindler 
regions, the continuation across from t = CT to t = + occurs via paths which run up (or 
down) the imaginary t axis and back again. On these paths, y is also evolving from — oo to 
+oo. Modes with nonzero momentum k y in the Milne region undergo an infinite number of 
oscillations as they approach t = 0~ from above or below, and an infinite number more as 
they cross the Rindler wedges. More subtle is the behaviour of the k y = modes. As we 
discussed, these modes generically diverge logarithmically as one approaches t — 0~. By a 
choice of phase one can put this divergence into the imaginary part of the mode functions. 
Then, as one circumnavigates the origin in the complex t-plane, the logarithm aquires an 
imaginary part of This causes the real part of the mode functions to undergo a jump, 
of just the amount needed to reverse its sign. This is illustrated in Figure 4. 
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FIG. 4. Behaviour of the k y = modes as they cross t = 0. With our choice of phase, the 
imaginary part (Im) diverges logarithmically but the real part (Re) is finite. Analytic continuation 
along the path shown in Figure 3 causes the real part to be odd in t whereas the imaginary part 
is even. 

The method described above is, we believe, completely adequate for dealing with quan- 
tum fields on A4 C - However, it is also interesting an important to develop the corresponding 
description of passage through the singularity in the d dimensional effective theory. In this 
theory, a scalar field has action 



with a specific time dependence in the kinetic term. Our approach here will be to regularize 
the theory by changing \t\ to It] 1 " 6 , with e a parameter analogous to that in dimensional 
regularization, to be taken to zero after renormalization. It is then necessary to add coun- 
terterms to the Hamiltonian at t — in order to render the time evolution operator well 
denned in the e — * limit. These counterterms have the effect of inducing a shift in the 
scalar field, proportional to its momentum, and analogous to the jump produced in the ana- 
lytic continuation method illustrated in Figure 4. This shift cancels divergences and renders 
the final state well defined. We shall show that within this method, demanding that the 
counterterms be local in x, and imposing time reversal symmetry is enough to uniquely fix 
the vacuum state. 

Finally, we shall point out an intriguing mapping between this problem and that of free 
fields on de Sitter spacetime. Under this mapping, the surface t = + corresponds to the past 
timelike infinity in de Sitter space, and t = CT corresponds to future timelike infinity. While 
these two surfaces are only connected at a point in the Milne universe, they are connected 
by a smooth bulk (comprising the entire spacetime) in the de Sitter case. The matching in 
de Sitter spacetime is unambiguous and again we shall show it corresponds to the previously 
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obtained results. There are holographic elements of this correspondence. Holography is 
naturally framed in terms of null surfaces,!^ and our approach involves matching information 
located on the two null surfaces t = 0~ and t = + . However, when we map to de Sitter 
spacetime, these two surfaces map to two spacelike surfaces, future and past timelike infinity, 
which are those which have been employed in the proposed de Sitter-CFT correspondence." 

All of these methods yield the same result for the quantum vacuum state on Aic- Be- 
cause there is no mixing of the positive and negative frequency modes, there is no particle 
production in the free field theory. However, once interactions are included, particle pro- 
duction occurs, and in Section VI we demonstrate that it is well defined. The k y ^ modes 
are produced with a density which tends to zero exponentially as H 5 vanishes, suggesting 
an adiabatic limit in which the particle production vanishes in the limit of slowly colliding 
boundary branes. The k y = modes do not show this behaviour, but we shall discuss how 
within string theory we can anticipate how an adiabatic limit may in fact emerge. 

Finally, let us mention the connection between this work and other, more ambitious at- 



tempts to direct 



y construct string theory on the compactified Milne spacetime considered 



here. Nekrasov^ considered string theory on the Lorentzian orbifold constructed by orbifold- 
ing Minkowski spacetime by a boost. In that construction, the two Rindler wedges become 
compactified in a timelike direction, producing two extra cones projecting horizontally from 
the origin, which possess closed timelike curves. Additionally, line segments emanating from 
the origin are produced in each of the four null directions. Cornalba and Costa avoid these 
features by modding out by a boost combined with a translation, replacing them instead 
by a new region containing a naked timelike singularity .0 Balasubramanian et al. consider 
other examples of time-dependent orbifold backgrounds in string theory. Whilst free strings 
seem to be well defined in these backgrounds, it is not yet clear whether interactions can be 
consistently introduced. 

The approach we suggest here does not amount to orbifolding Minkowski space. Instead, 
we use free field evolution (or, equivalently analytic continuation in t) to define a matching 
rule between the big crunch and the big bang. This difference is unimportant in the free 
theory, since the only difference in that case between our approach and the orbifold approach 
is that we would declare that the extra regions in the orbifold approaches do not exist. It 
is when we introduce interactions that the difference becomes crucial. In our case, the 
interaction Lagrangian is only integrated over the physical compactified Milne spacetime, 
whereas in the orbifold approaches it would be integrated over the additional regions too, 
containing closed time-like curves or naked singularities. We should stress that we have 
not attempted to construct string theory in our approach, therefore we cannot say whether 
string theory will ultimately be consistent on compactified Milne. However the field theory 
results are suggestive and we hope they will be a guide to such a construction. 

V. EMBEDDING MILNE IN MINKOWSKI 

A. Positive and Negative Frequency modes 

In this section we describe our first construction of quantum field theory on Aic x R d ^ x ■ 
A Fourier mode of a massless field, (p(t,y,x) = Lp(t)e^ k ' x+kyy \ obeys the field equation 
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where dot denotes partial derivative with respect to t. We have introduced the effective two 
dimensional mass m 2 = k 2 , and henceforth the x dependence shall play a purely spectator 
role. Equation (J7|) is just Bessel's equation, with imaginary order v = ik y /H 5 . It has a 
singular point at t = 0. The solutions which tend to positive and negative frequency WKB 
modes at late times are the Hankel functions, and the properly normalised outgoing positive 
and negative frequency modes are 



K = ^^^/ fls M^. r ky = ^f^H^ y/H5 {mt)e^y. (8) 

We would like to continue these modes to negative times. The Hankel functions have the 
following integral representations 



H$>{z) = e -4- 



which is analytic in the upper half z-plane, and 



/oo 
du exp(iz coshw — uu), (9) 



Hl 2 \z) = / duexp(— iz coshw — uu), (10) 

TTl J-oo 

which is analytic in the lower half z-plane. Consequently the Milne mode functions can be 
expressed as 



i / I f°° k 

^ = 2V ^j^dU^Vt-i™ 1 COShU ~ ijf^ ~ H M))' ^ky = ^ky- (11) 

By shifting the integration variable u — > u + yH^, we obtain 



i / I f°° k 
ipk = -\ — 77- / du exp( — iffit cosh u cosh H^y — ivrtt sinh u sinh H§y — i-^-u) . (12) 



(5 J-OO ±±5 

Further changing variables to Ky = — msinh(w) gives 



i I r°° dK y 



^ = 9\/^r/ I 7 Y , M Shlh ~ 1{KYM e*p(iK Y Y - ijKyT^T), (13) 
I V nH 5 J-00 kI + m 2 



where T = tcosh(H 5 y) and Y = tsmh(H 5 y) are the embedding co-ordinates in Minkowski 
space. This is a superposition of positive frequency plane wave modes on Minkowski space 
with momentum Ky. We note that the right hand side is an oscillatory integral which 
can be defined for t < by inserting a suitable convergence factor. Therefore the integral 
representation may be used as the definition of the mode functions there. 

The above integral representations of the Hankel functions define a natural analytic 
continuation across t = 0. One can read off from (§) and (|10"D the relations 

HW(e-™z) = -e^HjPiz), H^\e in z) = -e~ mv Hf\z). (14) 
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To see what these imply for the k y = modes, recall that 



H^\mt) = J (t) + zJV o (0, #o l V £ ) ee J (t) - jjV (t). (15) 

The rule (fL4l) implies that the analytic continuation of H^\mt) to negative values is 
—Jo(—t) — iNo(-t). Therefore from (|TJD the real part of Hq (mt) is an odd function of 
t, with a discontinuity at t = 0, and the imaginary part is even, with a logarithmic diver- 
gence at t — 0. The real and imaginary parts are illustrated in Figure 4. 

From the integral representations (§) and (|T0|) one can determine the behaviour of the 
analytically continued Hankel functions at large positive or negative t by performing the 
integral via the stationary phase method, obtaining 

*-±». (16) 

and similarly for Hj^ > H6 (mt), with i — > — i and fc^ — > — fcy. This continuation implies that 
there is no particle production since positive frequency incoming modes are matched to only 
positive frequency outgoing modes. 

It is important to stress that this choice of vacuum is priviledged. As we explained 
earlier if we cut off the singularity in a crude way, then for a generic choice of \in > and 
\out > states, for each k we would obtain in the k y = mode a particle production rate 
that diverges logarithmically with the regulator. It is only for the special case in which we 
define [out >— \in > or at least some finite Bogoliubov transformation of the \in > state 
that we obtain a finite result. From the d dimensional perspective this seems contrived, but 
from the d + 1 dimensional picture and the embedding in Minkowski space it is clearly the 
most natural choice. In later sections we shall also justify this matching from a purely d 
dimensional point of view. 

Finally, let us mention that our definition of in and out vacuum modes is not the same as 
that which has conventionally been used in treatments of quantum fields on Milne spacetime 
(see for example, Ref. pL8|) . In previous work, only the t > part of M. was used, and the 
initial vacuum was taken to be the 'conformal vacuum' as t — >• + , defined by the 'positive 
frequency' modes behaving as e~ lkyT as the conformal time r — > — oo. This is of course, not 
an adiabatic vacuum state, and therefore a somewhat arbitrary choice. In the conformal 
vacuum state, one finds particle production occurs in passing from the big bang t — > + to 
the asymptotic future, even in free field theory.0 

B. Projection onto M.c 

We have not yet distinguished between the Milne space M. and its compactification 
M.C which as we described above equation (^) is just A4/T(H 5 L), with T(H 5 L) the group 
of boosts with rapidity H 5 L. 

In Minkowski spacetime a particle is defined in a group theoretic sense as an irreducible 
projective representation of the Poincare group. We can similarly define particles on A4c by 
using representations on the covering Minkowski space that are invariant under the action 
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of the boost T{H$L). The map from Milne to Minkowski introduced in the previous section 
is inverted by means of a Fourier transform, to obtain 



* Ky (Y,T) = e^rr-i^T^T = j°° ^ e ^ H m /H ^ mt yk v {y- H ^^-\^)) (17) 

The plane waves ^K Y iXi T) form a representation of the two dimensional Poincare group, 
and the action of the boost A on these modes can be expressed as 

A : ^ Ky (Y, T) - V Ky (AY,AT) (18) 

where (AY, AT) = (tsmh(H5(y + L)),tcosh.(H^(y + L))) is simply a translation by Liny. A 
representation of the group Poincare /Y{H^L) can be constructed by simply summing over 
all boosts, 

oo oo 

q Ky (Y,T)= ^K Y (A n (H 5 L)Y,A n (H 5 L)T) = ]T * KY {A{nH 5 L)Y,A{nH 5 L)T). (19) 

n=— oo n=— oo 

We shall only use these functions on the physical region of interest, namely M.c, where they 
are given by 

* Ky (Y,T)= £ /" ^e^flg (^e*.^-^ 1 -^- 1 *). (20) 

Now using the Poisson summation formula J2^=-oo I^L dte 2mnt f(t) = J2m=-oo f ( m ) > we 
obtain 

^ y (Y,T) = i £ e^H^,H^y^^^-^\ (21) 

n=— oo 

This is just the expected result that summing over boosts projects out only those states that 
are translation invariant under y —>■ y + L, and is equivalent to quantizing the momentum 
k y = 27tn/L. If we were to perform the further projection onto the orbifold SyjZi mentioned 
in the introduction, we would now consider separating the k y modes into those which are 
odd and even under y — ► L — y. In string theory, this step introduces new states ('twisted 
states') but for field theory describing quantum mechanical particles, it has no such effect. 

The Feynman propagator on A4 is obtained by simply restricting the d + 1 dimensional 
Minkowski space propagator to the Milne region. In d + 1 dimensions the Feynman propa- 
gator iJH 

d-l 

G F (x,x>)=-i-^ ^) 4 Hfl iy2 (m{-a-ie)y% (22) 



where x = (T, Y, x) and a = (x — x' 2 ) + (Y — Y') 2 — (T — T') 2 . Restricting the Feynman 
propagator to the Milne patch simply requires writing a in terms of the Milne co-ordinates 
a(x, x') = (x — x 1 ) 2 — (t — t') 2 — Att' smh 2 (H 5 (y — y')/2). The Feynman propagator on the 
compactified Milne spacetime is obtained by projecting onto the boost invariant states. This 
is given by 
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oo 

G^ c (x,x') = £ GP(A n (H 5 L)x,x') 

n=—oo 

d-1 

= t -ijT^i^-) 4 H^_ 1)/2 {m{-a n -ie)^) (23) 
n=-oo (Am) 2 \a n + ie) 

where a n (x, x') = (x—x') 2 — (t—t') 2 —4tt'smh 2 (H 5 (y+nL—y')/2). The Feynman propagator, 
in addition to the interaction vertices, is all one needs in order compute the S matrix via 
perturbation theory on M.c- 

C. UV divergence behaviour 

It is important to understand whether compactifying Minkowski spacetime into M. c x 
R d ~ l introduces any new ultraviolet divergences. For the construction given above, the free 
field propagator on M. is just the Minkowski space propagator evaluated on M. Therefore 
it has just the usual divergences. In this section we shall show that the same is true for 
the propagator on Aic, for all points x and x' away from t — 0. This is to be expected 
intuitively since the Green functions on M. c are constructed by summing over boosts on 
one argument x, and these boosts carry x further and further from x'. 

The difference between the Feynman propagator on A4 C and M. is given in d + 1 dimen- 
sions by 

d-1 

oo / 2 \ ~4~ 

AG F (x,x') = £ -i—iTT HH H$ (m(-a n -iey/% (24) 
„=-oo,^o (4tu)— \<Tn + ieJ 

which in the coincidence limit x' — x becomes 

d-i 

oo / 2 \ ~4~ 

2 Y -i -itt o 27^ 7-, H$ 1 w 2 (m(4t 2 sinh 2 (# 5 nL/2) - le) 1 ' 2 ). 

ti (4ttz)^ \-4t 2 smh 2 (H 5 nL/2)+ie) (d " 1)/21 1 ^ 1 ' ' ' 

(25) 

The large \z\ asymptotic behaviour of the Hankel function is H^\z) — > (— ) 1 / 2 e - l ( 2 -3 I/7r_ i 7r ) 
and so the sum is rapidly convergent for nonzero t. Thus, at least away from t — the UV 
divergence behaviour of the Green function on M. c is just the same as that on Minkowski 
space. The behaviour of the Green function at t = is a more delicate matter, linked to the 
way in which interactions enter, which we shall discuss below. 

In the next section we shall see that if interactions are introduced on M. as integrals 
over fields on Ai, there are physical processes such as particle creation from the vacuum 
that occur at tree level, and which have no counterpart in Minkowski spacetime. They arise 
because energy is no longer conserved when the interactions are time- dependent. 

VI. INTERACTING FIELD THEORY 

The prescription discussed above for matching the big crunch phase to the big bang phase 
in the Milne universe relied on free field theory. That is, in the Minkowski spacetime within 
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which the Milne universe is embedded, we are propagating the fields according to the free 
field equations from the past light cone (on which Milne time is t = 0~) to the future light 
cone (on which Milne time is t — + ). With this prescription, as we have emphasized, there 
is no particle production. However, once interactions are included, particles are generically 
produced because the interaction terms in the Hamiltonian are time dependent. We shall 
calculate this effect in this section, using a very simple toy model of the interactions. This 
is not intended to accurately represent the actual interactions in string theory, but we hope 
will illustrate the general behaviour including the sensitivity to infrared and ultraviolet 
cutoffs. The former should come from cosmological evolution, since the growth of the extra 
dimension ceases when the universe becomes radiation dominated. Ultraviolet divergences 
have to be controlled by string theory or M theory effects and we shall comment on the 
possible form of these below. It is important to stress that we use the Minkowski embedding 
only to determine a matching condition in the free field theory. The interacting theory lives 
on the physical spacetime Aic- This is the sense in which our approach differs from one 
employing an interacting theory on a Lorentzian orbifold which is Minkowski space modulo 
a boost. 

As a very simple example, consider an interaction of the form S — — j ^/i 2 ip 2 , where 
the integral runs only over Aic- F° r concreteness we shall take d — 4, so there are three 
noncompact dimensions x. The interaction is simply a mass term, which from the point of 
view of the embedding theory in Minkowski space, is turned off outside the future and past 
light cones. We would like to compute the particle production due to this interaction. The 
quantum field ip is expanded in terms of creation and annihilation operators as 



where the k dependence of the modes is not explicitly shown. The creation and annihilation 
operators are normalized to obey 



a k y ,k> a i> y ,k> 



5 kv , K (27r) s 8 s (k-k'). (27) 



We may now compute the transition amplitude between the incoming vacuum state and 
an outgoing state with two particles, with equal and opposite momenta k y and k. The 
calculation is straightforwardly performed by first integrating over y to obtain the delta 
function 5k y +k y ,o corresponding to momentum conservation. Then we use the representation 
of the Hankel functions given in (|T3|) , evaluated at Y — 0, T = t, to obtain the interaction 
matrix element 



8Ha7rV 



2 

% -^—(2Tif5\k + k') J dt\t\ J J dudu'e^-^e-^^+^'^e-^, (28) 



where we have inserted a Lorentz invariant convergence factor, and include the volume factor 
needed to normalise the final states (since (2-7r) 3 £ 3 (0) = V). The integrals are straightfor- 
wardly performed using the identity 

00 , cosk y x Ak v it 

d x i = v. (29) 

-oo cosh 2 (a;/2) sinh(/cj / 7r/f/" 5 ) ' 
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to give 



<*>■ s - *• out| °- in> = 2^k? (a)3 ^<* + *> w (30) 

The probability for a transition from the vacuum to two particles with momenta k y = 
±27m/L and k within d 3 k is therefore 

£y ^k 1 (ky/Hsf 

16 (2tt) 2 |^|4sinh 2 (^7r/iJ 5 )' 

Dividing by the volume V one obtains the probability per unit volume for creating such 
particle pairs. At fixed external momentum, the final density of pairs is finite, as claimed. 
Furthermore, the k y 7^ modes which naively might be thought to be the most dangerous, 
are strongly suppressed. As the rate H 5 of contraction of the extra dimension is decreased, 
the production of these Kaluza Klein modes becomes exponentially small, showing the ex- 
istence of an adiabatic limit. The k y = modes do not however display such a limit, and 
in fact the result for particle creation for these modes is completely independent of H 5 . We 
shall discuss how this behaviour is likely to be altered in string theory, below. 

The integration over k in fl3T|) is infrared divergent. This is however an artefact of the 
fact that the interaction term we introduced diverges as t — > ±00. In the situations of 
interest for the cyclic and ekpyrotic models, the extra dimension tends to a maximum size 
t « and t » and this would introduce an infrared cutoff in \k\ of order Tq 1 where 
Tc is the characteristic time scale over which Milne-like behaviour holds. The total number 
density of created particles with this simple (p 2 interaction is ultraviolet finite. But the 
total energy density is logarithmically divergent. This disease may be cured by introducing 
the dilaton into the nonlinear field interactions, which generically occurs in string theory. 
We have e^v^ -2 )/^ -1 ) oc \t\, and each extra power of \t\ in the interaction introduces an 
extra negative power of \k\ in the matrix element, or \k\~ 2 in the probability. Conversely, 
if we introduce higher powers of (p in the interaction, e.g. (p 3 , this would boost the rate 
in the ultraviolet, just because more particles are created in each process and the phase 
space integral would involve a higher overall power of k. (Recall, there is no conservation 
of energy here since we have explicitly broken time translation invariance). Again, we can 
make the produced number or energy density finite by introducing sufficient powers of the 
dilaton coupling. However, as we shall now explain, we believe there should be additional 
effects suppressing the rate of production of particles with high momentum in string theory. 

The origin of the power law fall-off in these calculations may be traced to the sharpness 
in time t of the simple field theory interaction term we have introduced, and the fact that 
this interaction cannot be correct for very high momenta of the external particles. The 
Fourier transform of a function with a sharp kink falls off only as a power law, so with this 
field theory interaction, high energy energy-nonconserving processes are only suppressed as a 
power law of the energy. If interactions can be consistently introduced in string or M theory 
on A4c, we believe they will not show this sharp behaviour. Strings or membranes are never 
localised below a minimal length scale l s , and one would expect them not to see the extra 
dimension shrink below a minimal length scale l s . One could model this by replacing the 
factor H$\t\L, which is the length of the extra dimension, with Jl 2 + H 2 L 2 t 2 , which never 
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falls below l s . If we make this replacement in the particle production just computed, then 
the final particle density is actually exponentially convergent in the ultraviolet. Returning 
to (p8[), we see it is dominated by u ~ u' ~ 0. The effect of introducing the cutoff / s is 
therefore roughly the replacement 

j dt\t\e- imt ~ eltl -> J dtVt 2 + 5 2 e- imt - eltl , (32) 

where 5 = l s /(H 5 L). The left hand side equals — 2/m 2 as e — > 0. This exhibits the power 
law dependence of our answer for the amplitude above. However, for large mS the right 
hand side decays exponentially in m5. To see this, compute the difference between the left 
and right hand sides of fl32l), in which e may be set to zero from the outset. By integrating 
by parts, one can reduce the difference to 

2 f°° dt , , 5 2 

o+ / — -cos(mt) 3-. (33) 

m 2 Jo m 2 v l (t 2 + 5 2 )^ 

It follows that the integral equals minus the right hand side of (|32|). The latter is a Hankel 
function of imaginary argument, which decays exponentially, as e~ mS = e -mls/(HsL) ^ f Qr 
large ml s (i.e. particle momenta well above the string scale), or for small H 5 L (i.e. a small 
contraction speed of the extra dimension). 

Two important things occur in this model. First, the cutoff is not at the string scale, 
it is at |A;| = m ~ (H 5 L)/l s . More importantly, for fixed \k\, the particle production 
becomes exponentially small as H$L is lowered below \k\l s . This means that for H^L << 1 
and for modes of any fixed physical wavelength in the non-compact direction, there is an 
adiabatic limit in which the extra dimension can disappear and reappear with vanishingly 
small particle production. It remains to be seen whether these two desirable features will 
survive in a complete string theoretic calculation. 



VII. THE D DIMENSIONAL PERSPECTIVE 

The cyclic and ekpyrotic universe scenarios represent attempts at consistent cosmologies 
based on M and string theory. The resolution of what happens at big bang/big crunch 
singularities must be found within those theories. If, for example, the extra dimension 
involved is the eleventh dimension of M theory, as represented in the model of Horava- Witten, 
then when the two boundary branes approach the theory reduces to weakly coupled heterotic 
string theory. The full dynamics of the eleven dimensional theory involves eleven dimensional 
supergravity and the associated super-membrane. Nevertheless, one hopes that at least for 
slow motions, the system evolves quasi-adiabatically and at small brane separations one is in 
the regime of ten dimensional string theory. Furthermore, since the string coupling constant 
vanishes as the extra dimension disappears, string interactions should be suppressed at the 
singularity itself. However, as noted in Ref. the description of the bounce in string theory 
may not be straightforward, because the ten dimensional string frame metric vanishes at 
the boundary brane collision. 

In this section we shall attempt to describe not string theory but quantum field theory 
from a purely d dimensional perspective. As in the string theory case, the d dimensional 
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metric will vanish at the singularity. Nevertheless we shall see that the singularity may be 
traversed in a reasonably natural manner. 

The d dimensional Einstein frame geometry corresponding to the d + 1 Milne geometry 
considered earlier is given by 

ds 2 = \H 5 t\ 2 ^ d ~ 2 \-dt 2 + dx 2 ) (34) 

Note that the proper time t of the d + 1 geometry is identical to the conformal time t of the 
d dimensional Einstein frame geometry. The equation of motion of a massless scalar fieldP] 
is 

t 2 Cp +t(p + t 2 k 2 cp = 0, (35) 

just the equation for the k y = mode of a d + 1 dimensional massless scalar. Using the 
higher dimensional point of view as our guide we are lead to believe that the natural 'in' 
and 'out' states are 



r ut ' + (t,x) = ^ + (t,x) = ^— fr o ro(|&|)e** (36) 
and the corresponding Wightman function is 

G+(x,x') = J ^^^(|fc|(t-.e))^ 1) (|fc|(t + . e ))e^') (37) 

How can we understand why this is natural purely from a d dimensional point of view? As 
already discussed in Section III, regulating the spacetime in the sense of \t\ — > \/t 2 + 5 2 
will generically produce a particle production that diverges as hi(H 5 S). The fact that this 
divergence is identical for each k suggests that it may be removed by a counterterm which 
is local in d dimensions. 



A. What renormalization? 

How should we remove the logarithmic divergence in the scalar field as it approaches 
the singularity? From the d dimensional perspective, we would like to follow the traditional 
renormalization program: regularizing the theory, adding counterterms and finally removing 
the regulator. But we need to discuss what form the renormalization should take. 

As mentioned above, in the limit when t tends to zero, the dynamics of the field are 
dominated by the kinetic term. This term possesses a symmetry t — > Xt (which is just 
a translation in conformal time). The scalar field tends to ip ~ A + Bln\H^t\, and its 
momentum 7r tends to H 5 B. Rescaling time as above has the effect of increasing ip by BlnX, 
or i7^ 1 7rlnA. Therefore it is very natural to seek to exploit such a shift in order to remove 



Near t = massive particles behave like massless ones so it is sufficient to consider the massless 
case. 
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the divergence in ip. Since tt is asymptotically a constant, one can simply match it across 
t = 0. But to make (p finite, we need to redefine it via ip ^ p + Ctt, across the t = surface 
in the regulated theory, with the constant C chosen to obtain finite correlators for tp and tt 
for t > 0, as the regulator is removed. 

The shift <p —* <p + Ctt, tt —>■ tt, is a canonical transformation. It can be implemented 
by adding a a local counterterm to the Hamiltonian, which acts at t — to produce an 
additional unitary transformation taking the incoming to the outgoing quantum state. We 
shall show below that demanding the vacuum be both Hadamard and time reversal invariant 
uniquely fixes the value of the constant C. 

B. Dimensional Regular izat ion 

Before we can renormalize the <i-dimensional theory we must first regularize it. We shall 
use a regularization which makes both the field and its canonical momentum finite at t — 0, 
but which allows the background to remain a solution of the field equations everywhere 
except at t = 0. 

The dimensional regularization we use relies upon the generalization of the 1 + 1 dimen- 
sional Milne universe we have so far studied to a 1 + n dimensional Milne universe. Like 
their 1 + 1 dimensional cousin these are just a re- writing of Minkowski spacetime, but in 
dimensions n > 1 their constant time slices are negatively curved hyperboloids H n . We 
therefore consider the following vacuum Einstein's equations in d + 1 dimensions, 

d—n 

ds 2 = -dt 2 + H 2 t 2 dHl + Y, dx i, ( 38 ) 

where dH 2 is the line element on H„. This is a solution of the d + 1 dimensional field 
equations if the curvature (Ricci) scalar of the H n is —n(n — l)H 2 /2. From the d — n- 
dimensional point of view, we are considering fields which are constant on the H n , we have 
Einstein gravity plus a scalar field <fi representing the 'scale factor' of the H n .f\ The spatial 
curvature term in the d + 1 dimensional Einstein action leads to a potential for (f) in the 
dimensionally reduced d — n dimensional action, proportional to n(n — 1), and positive for 
n > 1. The idea of the regularization we use is to analytically continue this reduced theory 
in n, and ultimately take the limit as n — > 1. 

The action for a scalar field which is homogeneous on the H n is just 

- ~ J dt\H 5 t\ n V ^d^dv<p. (39) 



2 Parenthetically, we remark that even though H n is noncompact, there is a natural splitting 
between the homogeneous modes and the non-constant modes, because the Laplacian on H n has 
a smallest nonzero eigenvalue equal to the space curvature. Therefore provided we are considering 
d—n dimensional energy scales much smaller than H§, we may consistently neglect the non-constant 
modes. 



19 



For n = 1 — e, with e small and positive, it turns out that the field <p and its canonical 
momentum ir are both finite at t = and therefore both can be simultaneously matched 
across t = 0. To see this, note that as t — > the scalar field equation is approximated by 

t 2 Cp + (1 - e)t(p « 0, (40) 

with general solution (p ~ A + B(H^,t) e . The momentum conjugate to is given by 7r = 
|_f^ 5 t 1 1 ~ e <yi' ~ eH 5 B, constant in this limit. The remarkable feature is that for positive e, tp 
is also finite at t — 0, enabling us to match across £ = 0. (The limiting case e = 0, which 
we studied before, can be obtained by expanding (H 5 t) e ~ 1 + eln(if 5 £) for small e, and 
redefining the constants.) 

Whilst it is important that we have constructed the regularized backgrounds as solutions 
of the field equations everywhere except at £ = 0, for the purposes of this section all that 
really matters is that we have regularized the action for the scalar ip to (|39|) for t < and for 
t > so that we can match <p and 7r across £ = 0, with the introduction of local counterterms 
at that point. 



C. Matching Modes across the Singularity 

The regularized field equation for the ^/-independent modes is 

t 2 <p + (l-e)t(p + m 2 t 2 (p = 0, (41) 

with m 2 = k 2 as before. The general solution is ip = (mt)( <E//2 - ) x(mt), with x a Bessel function 
of order e/2. 

As t tends to zero, solutions to the regularized equation tend to the form ip ~ A + B\kt\ e , 
with A and i? constants. Thus as claimed above, both ip and 7r = iH^tl 1 ^ 6 ^ tend to constants 
as t tends to zero. Let (A~ , Z?~) and (A + , B + ) denote the values of these constants for t < 
and t > 0. Matching and 7r at £ = gives A + = A~ and i? + = —B~. Thus the asymptotic 
solutions for small \kt\ are 

V? ~ A + Bsign(t)\kt\ e . (42) 

Using the asymptotic form of the Hankel functions for small argument and redefining the 
constants we find that the general solution for all t is 

(p = 

\kt\^ 2 [C (e^ 2 H^ 2 (\kt\) + e-^ 2 H^{\kt\)) + Dsign(t) (H^(\kt\) + H^{\kt\))\. (43) 

As before defining the in and out states to be (pi n = A e / 2 H^l(—\k\t) for t < and (p out = 
A e /2H^2(\k\t) for t > we find the following relation between the in and out states: 

<p in = a<p out + (3<p* out , a = -icosec(ne/2)e~ lne/2 , (3 = -icot(ire/2). (44) 

Let us consider the limit e — > 0. In this limit a — ► — 2i/(ixe) and /? — > —2i/(ne). This 
implies infinite particle production, as we obtained before with the simple cutoff. In fact 
the analogy is 
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(45) 

This relationship is similar to that obtained in ordinary renormalization where l/(d — 4) 
divergences in dimensional regularization corresponds to the ln(A) divergences obtained with 
a UV cutoff. 



D. Regularization Independence 

In the previous section we have seen that the logarithmic divergences of the naive cutoff 
regularization used earlier manifest themselves as 1/e divergences in the dimensional regular- 
ization scheme. We now want to renormalize these divergences, by absorbing the divergent 
pieces in local counterterms added to the action. 

In dimensional regularization the minimal subtraction scheme is simply to throw away 
1/e divergences. Here, one must be more careful since there is a danger of violating unitarity. 
We are able to perform the subtraction via a unitary transformation as follows. First write 
the Bogoliubov transformation in the form 

-icosec(ire/2)e~ i7Te/2 -icot(vre/2) \ / <p in \ . , 
icot(7re/2) icosec(7re/2)e^ e/2 ) \ <p* n J ' 1 } 

Now redefine the out state by means of the following Bogoliubov transformation 

(47) 







Then in the limit e — > we find 



Vout 

/* 
'Pout 




(48) 



but since the 1/e divergences have been removed by means of a Bogoliubov (unitary) trans- 
formation, unitarity is clearly preserved. 

The overall minus sign aquired by the mode functions ip' out after passing through t — is 
an unobservable phase in the quantum mechanical wavefunction. However, it is interesting 
to compare what we have done here to what was done in the earlier methods of analytic 
continuation. Here we have insisted that 7r be continuous, and have found that, after renor- 
malization, (p must change sign. Whereas there, the continuation produced a divergent part 
of the field which was even with a corresponding momentum n which was odd, and had a 
jump at t — 0. The two methods differ by an overall minus sign but this is just a phase and 
is physically irrelevant. 

In fact, we can easily compute the magnitude of the jump in (p in the regulated theory, 
using the small argument form for the Hankel functions. We find the result that after the 
above renormalization, 

V?(0 + )-v?(0-) = -J-7r(0). (49) 
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But this is precisely a canonical transformation corresponding to the shift symmetry of the 
massless homogeneous scalar field described earlier. This justifies the conjecture than the 
divergent piece can be removed by means of a canonical tranformation. The result is that 
this subtraction scheme predicts no particle creation for a massless scalar field, exactly the 
result of the Minkowski embedding matching. So we have reached the same conclusion from 
the d-dimensional persepective. 

E. Unitary Evolution and Regulating the Hamiltonian 

The previous result is in a sense a trivial one. We can always redefine the in state 
by means of a Bogoliubov transformation such that [out >= \in >. However, the critical 
point is that the Bogoliubov transformation we used was momentum independent and so 
corresponds to a local unitary transformation. Any canonical/Bogoliubov transformation 
can be thought of as an instantaneous interaction term in the Hamiltonian. Let AS denote 
the 'impulse' (which has dimensions of action), acting at time t, so the time dependent 
Hamiltonian is 

H = H + 6(t)AS. (50) 

Then the unitary evolution operator Ujjit, t') for t' < and t > may be expressed formally 
as 

U H (t,t') =Texp(-iJ*H(t)dt) = U H (t, + ) exp(-iAS)U H (0~ , t'), (51) 

in other words unitary evolution with Hamiltonian Hq up until time t — 0~, a unitary 
transformation exp(— iAS) followed by unitary evolution with Hq again. In particular if AS* 
is quadratic in the fields and momenta then exp(— iAS) performs a canonical transformation 
of the form described in Appendix 1. Comparing with the result ( |49D above, we see that 
AS* is the local counterterm 

AS = — f d 3 x7r 2 (x,0). (52) 

eH 5 J 

In principle we have freedom to add other local counterterms, a freedom leaving us with a 3 
parameter family of matching conditions. However, as we explain in Appendix 2. the above 
choice is the unique one that yields a Green function which is both of Hadamard form and 
time reversal invariant, i.e. such that the Wightman function transforms under the time 
reversal operator T by G + (Tx,Tx') = G~(x,x'). 

Another way of seeing why AS" given in ( j52l) is the correct counterterm is to consider the 
bare Hamiltonian in the dimensional regulation scheme. Near the singularity the leading 
contribution to the Hamiltonian is from the kinetic term. Computing from the action, 
(P§|), we find the kinetic term contributes H = | / d 3 x7i 2 /\H 5 t\ 1 ~ t . The unitary evolution 
operator depends on time ordered products of integrals of this Hamiltonian. Consider the 
contribution of the region in the vicinity of t — 0, to / Hdt for \t\ < T. Using the fact that 
7r is nearly constant, we may approximate this for small T as 

- f T \H 5 t\~^dt I d 3 xn 2 {x,0) = iElDl f d 3 x7c 2 {x,0). (53) 
2 J~t J eHfy J 
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Now taking the limit as e — > 0, we obtain a leading term of precisely minus AS given in 
fl52f). Thus we see that the counterterm we have introduced is just the minimal subtraction 
needed to obtain a well defined limit for the time evolution operator as e is taken to zero. 

Let us briefly discuss the issue of general covariance. In this paper we are working in 
a fixed background and not including gravitational effects. There is a preferred slicing of 
this background and therefore it is allowable to have a counterterm which is explicitly t- 
dependent. However, when gravity is included, there should be no preferred coordinate. 
How should we think of the counterterm AS in that case? The point is that from the d 
dimensional point of view, there is a scalar field which yields a coordinate invariant time- 
slicing of the geometry, namely the dilaton <p which we have mostly neglected. If we rewrite 
AS" in terms of the dilaton field, it is then coordinate invariant and can be used in a generally 
covariant treatment. 



VIII. QUANTUM FIELDS IN DE SITTER SPACE 

In this section we want to discuss an interesting formal correspondence between quantum 
fields on M.c and on de Sitter space. The spacetime M.c x R d ~ l is conformal to the non- 
singular spacetime S 1 x dS d , as may immediately be seen from the metric 

ds 2 = Hlt 2 [dy 2 + — — ^ {-dt 2 + dx 2 )}, (54) 

[H 5 ty 

where the latter term is just the metric on de Sitter space in the flat slicing. The parameter 
H5 plays the role of the Hubble constant of the d dimensional de Sitter space. Figure 5 shows 
the correspondence between the two halves of the Milne universe, L and U, and the two 
flat-sliced regions of de Sitter space, one of which covers the region containing past timelike 
infinity, t = + up to the co-ordinate horizon at t = +00, and the other which covers the 
region from the horizon t = —00 to future timelike infinity t — 0~. In the de Sitter geometry, 
starting from a generic point in region L, it takes an infinite proper time to reach t = 0~, 
but in M. only a finite proper time is needed. The key point in this correspondence is that 
the two surfaces t = 0~ and t = + meeting at the singularity of M.c are mapped to the 
two spacelike surfaces representing timelike past and future infinity in de Sitter space. 
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FIG. 5. The Milne spacetime is both locally flat (left diagram) and conformally S 1 x de Sitter 
(right diagram, with the 5th dimension suppressed). In the de Sitter picture passing through the 
singularity corresponds to matching future timelike infinity in de Sitter to past timelike infinity. 

A peculiarity of this map is that the collapsing region L of the Milne geometry is mapped 
to the expanding region L of de Sitter space. Nevertheless the arrow of time always points 
in the direction of increasing t. 

De Sitter space can also be globally covered by the closed slicing co-ordinates, so-called 
because the Cauchy surfaces are spheres S' d_1 . If d 2 Q S d-i denotes the metric on an S^ 1 
with radius H$ then the de Sitter metric is 

ds 2 = -dr 2 + cosh 2 (H 5 T)d 2 n sd -i (55) 

There is no co-ordinate singularity in these co-ordinates which provide an unambiguous 
method for matching fields from Z~ to X + . There exists a unitary operator U{t\, t 2 ) which 
generates time evolution from r 2 to T\. This operator satisfies the Schrodinger equation 

i-^U(n,r 2 ) = H(r 1 )U(r 1 ,r 2 ) (56) 

where H(t) is the time-dependent Hamiltonian. The S matrix is denned as the unitary 
operator 

S= lim U(t x ,t 2 ) = Te- l ^ H(T)dT , (57) 

which in the Milne picture corresponds to an S matrix matching from t = + to t = 0~ or 
equivalently, 

|^,* = + >=S ,t ftM = (T> (58) 

Free field theory on de Sitter spacetime therefore provides us with yet another matching 
prescription for Milne spacetime. 
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Since the spacetime M.c x R d ~ l is locally flat then a massless minimally coupled scalar 
field is identical to a massless conformally coupled scalar field. In d + 1 dimensions the 
conformal coupling term in the Lagrangian is — Rip 2 . After a conformal transformation 
we obtain a massless conformally coupled scalar field on S 1 x dS d . The Ricci scalar on 
S 1 x dS d is simply a constant R = d(d — 1)H 2 which means we may reinterpret the scalar 
field as a minimally coupled scalar field with mass m 2 = (d — 1) 2 H 2 /A. On Kaluza-Klein 
compactification over the S 1 we obtain a tower of scalar fields on dS d with masses m 2 = 
(d-l) 2 H 2 /A + (2nn/L) 2 . 

As we explained earlier the KK zero mode of a minimally coupled scalar on M. c x R d ~ x 
is identical to a minimally coupled scalar on the lower dimensional Einstein frame geometry. 
This in turn is equivalent to a scalar field of mass m 2 = (d — 1) 2 H 2 /A on dS d . This result 
would seem surprising had we not used the higher dimensional geometry since although 
the d dimensional geometry is conformal to dS d , a minimally coupled scalar field is not 
conformally invariant. 

In de Sitter space there is a one-parameter family of de Sitter invariant vacua. E A 
unique choice, known as the Bunch-Davies vacuum,Ej is obtained if we also demand that 
the Feynman propagator be of the Hadamard form (see Appendix 2). A de Sitter invariant 
measure of the distance between two points is provided by the variable z = 1 + ((£ — t') 2 — 
(x — x') 2 )/2tt'. Although we have expressed it in terms of flat slicing co-ordinates, z is 
globally defined and has the same definition regardless of which patch t and i! are in. This 
means that the Feynman propagator in the Bunch-Davies vacuum can be defined globally 
in terms of the hypergeometric functional 

or((d-l)/2 + i/)r((d-l)/2-i/) , 1, 
G F = H 2 ^ » A-^ >1 i 2 F x [{d - l)/2 -u,(d- l)/2 + v;2,-(l+z- is)] 

(59) 

where v depends on the Kaluza-Klein mode number n, v = 2irin/ H§L. 

We are interested in calculating the S matrix to determine the matching condition implied 
by the above propagator. At the free field level it is sufficient to consider how a single particle 
state evolves. The Feynman propagator allows us to evolve the positive frequency part of a 
scalar field from one Cauchy surface t' = constant, to any point in its causal future, 



(x, t)=if d d - 1 x'(Ht , y d G F (x, t; x', t')($, - %)^> + {x\ t') (60) 



A simple calculation shows that a positive frequency WKB in state in the region t — > +oo 
evolves to a positive frequency WKB state in the region t — > — oo. The reason for this can 
be seen immediately from the conformal diagram since 1 + and X~ are identified as the same 
coordinate horizon. 

This result can be understood more clearly by realizing that in analogy with Minkowski, 
we can define particles in de Sitter as representations of the de Sitter group. Since these 
representations are globally defined and we choose a vacuum that respects the de Sitter 
symmetry, then it is clear that in the absence of interactions there can be no particle creation 
in de Sitter spacetime. 

At this point we should make clear that this is not in contradiction with the usual 
statement that there is a thermal distribution of particles in de Sitter space. This description 



25 



arises because of a different definition of particles, one which is more appropriate to the 
static patch surrounding an observer's world-line. The definition of particles in a curved 
spacetime is observer dependent, but the evolution of fields is observer-independent The 
Bunch-Davies vacuum is the vacuum in which globally there is no particle creation according 
to the representation theory definition, but where locally an observer sees a thermal bath of 
particles. 

In conclusion we have reached the same results as before for free field theory, using the de 
Sitter picture. When one includes interactions, a careful track of the non-minimal couplings 
must be taken into account. For instance \(p A theory in the 4d Einstein frame geometry 
will correspond to massive X'(t)(p 4 theory on de Sitter where X'(t) is a new time-dependent 
coupling constant. Ultimately, it is not clear to us that the Milne-de Sitter correspondence 
be a useful guide, because the problem of string theory on de Sitter space is probably a 
harder problem than that of string theory on the Milne spacetime. 



IX. CONCLUSION 

In this paper we have shown that it is possible to define free quantum fields on the 
compactified Milne universe in a consistent and unambiguous manner. We have made limited 
progress in studying interactions, and how these lead to particle production. The density of 
particles produced at fixed external momentum is finite at tree level. The integrated density 
was also found to be finite provided the dilaton dependence of couplings caused them to 
vanish sufficiently rapidly with t. We suggested how an adiabatic limit, in which particle 
production would be exponentially small for small H 5 L, might emerge in string theory. We 
also pointed out connections with quantum field theory on de Sitter spacetime, which may 
well be interesting in their own right. 

Certainly, much remains to be done to explore quantum fields on the Milne and com- 
pactified Milne universe. The methods used here could be extended to include gravitational 
backreaction, at least for linearised gravity, to follow cosmological perturbations through the 
singularity. We shall report on a study of loop diagrams for scalar field interactions on M. c 
in the near future. A major challenge remaining is to extend these ideas to string theory 
and M theory. 
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Note Added: Since this article was submitted to the archive a number of related articles 
condisidering fields/strings on backgrounds with cosmological singularities have appeared.0 
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X. APPENDIX 1: THE RELATIONSHIP BETWEEN CANONICAL AND 
BOGOLIUBOV TRANSFORMATIONS 



The group of Bogoliubov transformations is identical to the group of linear canonical 
transformations. For a 2 dimensional phase space these form the group Sp(2), ie those 2x2 
matrices that satisfy 

A T VLA = n, (61) 

where f2 = ^ ^ (/ ) ' * s a rea ^ ^ orm °^ (2) and consequently we can write any 

generator A in terms of the Pauli matrices a± = |(<7i ± cr 2 ), a 3 . 

A = A + {a)A_{b)A 3 {c) = e aa+ e ba - e ca * (62) 
If ^ ^ ^ denotes an arbitrary vector in phase space then the transformed vector ^ ^, ^ = 
A I ^ ] is given by 



' q \ = ( q 

p J yp + bq 
q \ _ ( e c q 



A. 



MpJ = U-< p J- <63) 

In quantum mechanics q and p are replaced by operators Q and P satisfying the Heisenberg 
algebra 

[Q,Q] = [P,P] = 0, [Q,P]=i. (64) 

The canonical transformations A + , A_ and A 3 can be represented by 3 unitary transforma- 
tions U+, U- and U 3 given by 

U+(a) = exp(-iaP 2 ), 
2 

U_(b) = e W (-UbQ 2 ), 



U 3 (c)=exp(^c(PQ + QP)). 



(65) 



These follow as a simple consequence of the Heisenberg algebra. We can define creation and 
annihilation operators in the usual way a = -j^{Q + iP), at = (Q — iP) and consequently 
an arbitrary canonical transformation corresponds to an arbitrary redefinition of a and at , 
ie. a Bogoliubov transformation. 
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In field theory an operator valued field can be expressed in terms of creation and anni- 
hilation operators as 



<p(x) = Y, a M x ) + Q k*(z) = V2Y^QiRe(ipi(x)) - Pilmtyiix)) (66) 

i i 

with ipi(x) and ip*(x) the positive and negative frequency mode functions, normalised ac- 
cording to the Klein-Gordon inner product. It is usual to write a Bogoliubov transformation 
as a transformation acting on these modes: ^[(x) = aipi(x) + (3ip*(x), which preserves the 
Klein-Gordon norm if |«| 2 — \(3\ 2 = 1. Since the field <p(x) is invariant under this transforma- 
tion, the creation and annihilation operators must transform under the inverse Bogoliubov 
transformation, 

a[ = a*cii — (3*a\ 

d- = aa\ — (67) 

If we re-write this in terms of coordinates and momenta, then in the above notation, the 
Boguliubov transformation corresponds to the linear canonical transformation on the oper- 
ators (Q u Pi) 

( Re{a - (3) Im(a + f3) \ . . 

\Im{{3-a) Re{a + (3) J 1 } 

It is simple to check that detA = \a\ 2 — \(3\ 2 = 1 as required. This formula gives the precise 
map between canonical transformations and a Bogoliubov transformations. 



XI. APPENDIX 2: HADAMARD FORM OF THE PROPAGATOR 

In the description of quantum fields on a curved spacetime, a natural and common 
restriction on the choice of vacuum is to impose that the Feynman propagator should be of 
Hadamard form. More precisely this means that in the coincidence limit x — > x', Gf(x,x') 
has the same singularity structure as the Feynman propagator on flat space. A physical 
motivation for this choice of vacuum is that two observers located at nearby points x and x' 
should not be able to tell if they are on a curved space or Minkowski space by information 
sent between them. Since we are considering the coincidence limit x — > x' the distinction 
between the various types of Green's functions is not relevant and it is common to work 
exclusively with the Hadamard function defined by 

G H (x, x') = G + (x, x') + G~(x, x') = -2lmG F (x, x'), (69) 

where G + (x, x') denotes the Wightman function < 0\{p(x)(p(x')\0 > and G~(x, x') is its 
complex conjugate. Now suppose we are interested in the Green's functions on a spacetime 
which is invariant under the discrete symmetries of time reversal T and parity P. This is true 
of all the spacetimes we have considered which can be seen by a simple inspection of their 
metrics. In particular then the spacetimes are invariant under the combined symmetries 
PT. This operator is anti-unitary and modes ipi(x) may be decomposed into its eigenstates, 
which with an appropriate choice of phase can be chosen to obey 
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ip(PTx) =i>*(x). 



(70) 



Given one such eigenstate, let us construct an arbitrary momentum-independent Bogoliubov 
transformation, to the mode function 

rj}>{x) = aip(x) +(3i/j*(x). (71) 

Then we find that 

iP'(PTx) = aip*(x) + [3ip(x). (72) 

If we demand that ip'(x) also be an eigenstate of PT with eigenvalue 77 (which can be an 
arbitrary complex phase), then we must have 

a* = i]a,P* = r/P, (73) 

relations which can only be satisfied if Im(af3*) is zero, a result we shall use in a moment. 

Given a particular set of positive frequency modes ipi(x) which are time reversal invariant 
ipi(PTx) = ipi(x)* and for which the vacuum is Hadamard we can construct the Hadamard 
function as, 

G H (x,x') = X>iO*0$V) + c - c - ( 74 ) 

i 

Now define a new vacuum by means of a constant Bogoliubov transformation of the first 
vacuum 

rt(x)=a1>i(x)+/3tf(x). (75) 

The Hadamard function in the new vacuum is given by 

G' H (x,x') = (|a| 2 + \(3\ 2 )G H (x,x') + 2Re(a(3*)G H (x,PTx') - 2Im(a(3*)A(x, PTx') 

The commutator function 

[<p(x), <p(x')] = G + (x, x') — G~(x, x') = iA(x, x') (76) 

does not contribute to the Hadamard singularity structure since it is vacuum independent. 
The singularity of A(x,x') occurs for null separated points only. However, the commutator 
G' H (x, x') in the new vacuum will have singular behaviour in the coincidence limit x' — > PTx. 
But one of the requirements of the Hadamard vacuum is that the propagator only has a 
singularity as x — > x' . Consequently if we demand that the new vacuum is Hadamard then 
we must have Re(a ( 5*) = 0. 

The conditions derived in the previous two paragraphs together imply that (3 — 0. In 
other words, the requirement of PT invariance and Hadamard form uniquely picks out the 
vacuum. In fact all the vacua we have considered in this paper are trivially P invariant, 
G + (Px, Px') = G + (x,x'), and so we only need to additionally impose the restriction of T 
invariance. 

To put this in the more familiar setting of quantum fields on de Sitter space, the condi- 
tion that the vacuum propagator be de Sitter invariant automatically picks out a vacuum 
invariant under PT since PT is a discrete subgroup of the full de Sitter symmetry. Then, 
as is well known, the additional requirement of Hadamard form uniquely picks the vacuum 
as the standard Euclidean, or Bunch-Davies vacuum. 
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